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In this paper, we explain the fundamental properties of the radiation 

processess of untrapped kinks moving on discrete lattices or any spatially 

periodic potential. In particular we explain qualitatively and quantitatively 

the interesting recent experiments of small discrete periodic dynamical arrays 

of Josephson junctions. We obtain and show the resonance condition for a 

kink to radiate phonons is the kink frequency or one of its harmonics KjX 

where X is the velocity of the kink and the Kj are vectors of the reciprocal 

lattice must equal one of the phonon frequencies, u)(k), where k are the 

lattice vectors determined by the lattice spacing and the periodic boundary 

conditions. Then we show the radiation resonance conditions determine the 

voltage at line center in the experimental currently voltage (J — V) diagrams. 

We also show the phonon linewidth due to external damping and radiation 

damping gives rise to steps and hysteresis in the I — V diagrams. In addition 

in the cases of more than one kink, say M kinks, there is a simple allowed 

value of k in contrast to the M — 1 case where there are at least as many 

resonances as there are phonon modes. We show the unique ku for a given 

M is the allowed k such that the distance between the kinks on the ring is 

one wavelength i.e. Xm = (27r)/Q . 



I. INTRODUCTION 

The purpose of this paper is to explain the fundamental properties of the 
radiation process of untrapped kinks moving on discrete lattices or any spa- 
cially, periodic potential. One of our principle aims is to explain the very 
interesting experiments and simulations of small discrete periodic dynamical 
arrays of Josephson junctions in a recent paper by Watanabe et al. [1]. In 
section II we first explain the radiation properties of kinks in the large (quasi- 
continuous phonon spectra) discrete sine-Gordon (SG) system because most 
of the fundamental kink radiation interactions are quite similar in the large 
and small systems. The solution of the kink radiation problem for the large 
discrete SG system by the rigorous collective variable (CV) theory including 
radiation from trapped kinks of [2] was found to be in excellent agreement 
with the simulations of the discrete SG even in the strongly nonperturbative 
regime. Unfortunately the full nonperturbative equations do not lead to a 
simple description of the radiative interactions. However when a perturbative 
treatment of the kink radiation process is valid we find we can understand 
the results in a simple intuitive way. The circular arrays of Josphson junction 
experiments are in a parameter range where the kink radiation interaction 
is treatable by perturbation theory between the center of the mass of the 
discrete SG kink X(t), which satisfies the non linear Peirels-Nabarro (PN) 
equation of motion in the presence of the kink. The radiation phenomena 
that occur in the Josephson junction experiments and in the large discrete 
SG systems we analyze in this paper come from kinks that are not trapped 
in the PN but are almost free i.e., the kink kinetic energy is much larger than 
the kink potential energy. 

An important early large discrete system radiation result for which per- 
turbation theory is valid was that of [3] who simulated and analyzed a large 
discrete SG system and found for a variety of initial conditions that there 
were certain special velocities of the kink at which the rate of loss of energy 
due to radiation of phonons changed abruptly. They called these special ve- 
locities "knees." They obtained numerical values for the velocity knees by 
a graphical solution of the dispersion law of the discrete SG which agreed 
with the simulations. In [2] we pointed out the velocity of the kink at the 
knees satisfied the condition (27ijX/a) = KjX = ujq, where uoq is the lower 
band edge of the phonon band, X is the velocity of the center of mass of the 
kink,a is the lattice spacing and j is an integer. In the present paper we show 
the radiation resonance conditions, KjX = ujq, for the large system is a spe- 



cial case of the general radiation resonance conditions which play the crucial 
role determining the value of the voltages in the current-voltage characteris- 
tics in the small discrete periodic systems. The frequency KjX is the j'-th 
harmonic of 2-nX/a, which we call the kink frequency, for reasons that will 
become clear in Sec. II. In the small discrete system there are a set of radi- 
ation resonance condition KjX /a = uj{k) where u(k) = [u>q + 4:sin 2 ka/2} 1 ^ 2 , 
k is the wave vector of the kth phonon mode, and the Kj are the reciprocal 
lattice vectors of our discrete lattice spacing of a. We show the large and 
small discrete systems radiation resonance conditions arise from the calcula- 
tion of the radiation from a kink moving at approximately constant velocity, 
X. In the large system the radiation resonance conditions give directly the 
value of the velocity knees Xj = oj§K~ x a. One of the conditions for the 
validity of perturbation theory of kink radiation is different in the large and 
small system cases. In the large system the phonon states form essentially a 
continuum of phonon states which leads to the requirement, uj t t ^> 1, where 
r r is the radiation relaxation-time. This condition is satisfied in [3] and all 
the large systems considered in this paper. While in the small systems where 
because of the large spacing between phonon frequencies it is necessary to 
have a weak external dissipation term, T0, in the discrete SG equation to 
provide the phonon line with a linewidth which leads to the requirements, 
Tr r ^> 1 which is satisfied in the experiments. In addition to T in the small 
system there is a contribution to the phonon linewidth due to spontaneous 
emission, [t(A;)] _1 , where r is the phonon lifetime which is a non linear func- 
tion of k and X. We gain further insight into the kink radiation problem 
by showing how the continuum SG with a spacially periodic perturbation, 
cos k'x 4>(x, t), which has been solved in perturbation theory and has a single 
knee at k'X = uq. This example makes clear the importance of spatial pe- 
riodicity of the potential in the kink radiation problem by defining the kink 
frequency, k'X in the perturbed continuum. SG and KjX in the discrete 
periodic lattice case. 

In Sec. Ill, we show the number of degrees of freedom and the periodic 
boundary conditions determine the allowed wave vectors which together with 
the radiation resonance conditions determine the velocities of the kinks at 
the phonon line centers. The voltage is proportional to X which in turn 
is a non linear function of the driver /. Consequently we find / is a non 
linear function of V which leads to the steps in the I — V diagrams. The 
non-linear dependence of I on V in small driven damped periodic rings leads 
to hysteresis between adjacent steps. 



We find the perturbative treatment of the CV equations of motion explain 
quantitatively the main experimental and numerical results of Watanabe et 
al. [1] on rings of Josephson junctions who consider rings of N = 4 and N = 8 
junctions with M — 1, 2, 3 and 4 kinks for various values of the discreteness 
parameter. In all cases the perturbation theory is valid and agrees with the 
experiments and simulations to within three percent. In addition in the cases 
of more than one kink, say M kinks, there is a simple allowed value of k in 
contrast to the M = 1 case where there are at least as many resonances as 
there are phonon modes. We show the unique /cm for a given M is the allowed 
k such that the distance between the kinks on the ring is one wavelength i.e. 
A M = (^)kM 1 . 

II. RADIATION-RESONANCE CONDITIONS IN LARGE 
AND SMALL SYSTEMS 

In this section we analyze the radiation mechanism of almost freely moving 
kinks in large and small discrete lattices where the perturbative effect of 
discreteness causes the kinks to radiate phonons. The continuum SG in 
dimensionless variables is : 

0-0" + (7r// o ) 2 sin^ = O (2.1) 

where the dot represents time derivative, the prime represents spatial 
derivative and (ir/l ) 2 is the dimensionless constant. The size of the kink in 
these variables ~ l . The dimensionless parameter ir/l is the slope of the 
kink and (uj = ir/l ) is the frequency of the lower phonon band edge. A 
complete discussion of the units is given in [2]. The discrete SG in the same 
units with the lattice spacing a = 1 is : 

<j> n - A 2 n + (tt// ) 2 sin 0„ = (2.2) 

where 

^2^-n = A n+ i + A n -i — 2A n 

The values (tt/1 ) > 1 are considered discrete and the continuum limit 
corresponds to n/l <^ 1. We assume a solution in the form <p n = f n + q n 
where f n = 4tan _1 exp[(7r// )(n — X(t))] and Ti 2 qjfj = and J2Pjfj = 
are the constraints where /' is the derivative of /_,- with respect to X(t) and 
is called the shape mode. The derivation of the complete CV equations of 



motion for q n (t) and X(t) are given in [2.] In all cases we consider the kink 
radiation interaction can be treated perturbatively. Consequently we need 
only the lowest order theory namely the nonlinear equation for the center of 
the mass X(t) in the Peirels-Nabarro PN potential and the linear equations 
for the phonons q n on the presence of the kink. The lowest order equation of 
motion for X(t) is 

M(X)X = -J2 £[A 2 / B - {-n/lof sin f n ] + \x 2 ^- Y.U'nf 

and terms linear in q see Eq. (2.11) of [2], where M(x) = J2{fn) 2 = M + 
M(X) = 8(tt// ) + AM(X). 

For the levels of discreteness in the experiments and simulations we ana- 
lyze in this paper AM(X) <C M thus we can replace M(X) by M and we 
can neglect the X 2 term. The linear terms in q in the X equation represent 
induced absorption and emission terms which are negligible compared with 
the spontaneous emission term in the experiments and simulations analyzed 
in this paper. Finally we obtain 



X = -{M )- l Y,f n [^2fn ~ (7r/g 2 sin/„] = ]T C j sin 2*jX(t)/a (2.3) 

n 3 

The dominant term in Cj goes as Cj exp(— jirl ). Thus Cj decreases with 
increasing j and l . Thus as the size of the kink, 2/ , increases compared with 
a, Cj — > and in this limit of the continuum the force due to discreteness 
vanishes. For appreciable discreteness effects the size of the kink should be 
less then six or seven lattice spacings. The linearized equations of motion for 
the q n are: 

q n - A 2 q n + [(vr// ) 2 cos 2irf n ]q n = Xf' n + A 2 /„ - (vr// ) 2 sin 2nf n 

= ^^(n)sin27rjX/a (2.4) 

3 

The closed coupled set of equations for X and q n Eqs. (2.3) and (2.4) 
constitute the lowest order CV equations which explain the spontaneous 
emission of a SG kink in a discrete lattice. The normal modes and eigenvalues 



for the q n when the right hand side of Eqs. (2.3) and (2.4) are set equal to 
zero are 

uj 2 (k) = (n/l ) 2 + 4sm 2 (k/2) 

and 

il>k(n) = [27ra;(/c)] _1 exp[A;7rn// ]{iA; — tanh[7m/7 ]} (2-5) 

The phonon modes in the presence of the kink, Eq. (2.5) are orthogonal to 
the shape mode of the kink f' n . Note the spectra u(k) is the same whether 
the kink is present or not. Strictly speaking Eq. (2.5) and u(k) are rig- 
orously true in the continuum limit and approximately true in the weakly 
and moderately discrete region. All of the experiments and simulations we 
consider in this paper are in the moderately and weakly discrete limit. In the 
extreme discrete limit n/l > 1 where the kink is only two or three lattice 
spacings the eigenfunction and eigenvalues are not known analytically but 
only numerically because in the extreme discrete limit the ground state is 
only known numerically. We thus have a model of a point "particle" X(t) 
with Hamiltonian 

H - = ik + Vf " (x) (2 - 6) 

where 

V PN {X) = 2Ci(2^') _1 [l " cos2irjX/a] 

spontaneously emitting phonons through the terms on the right hand 
side of Eq. (2.4). The PN frequency up N = J2C 2 2ttj. From perturbation 
theory the energy of the kink plus the energy of the phonons is conserved. 

There are two kinds of spontaneous emission of phonons. The first comes 

p 2 
from a kink trapped in a PN well, i.e., when j^- < Vpn(X) and the kink 

radiates harmonics of ujpn where mupN > ou and m is an integer. This is 

anharmonic radiation which is treated in detail in [2] and plays no role in 

the phenomena analyzed in this paper. The phonon radiation mechanism 

which explains the phenomena analyzed in this paper is determined by the 

frequency which the center of mass of the kink, X(t), passes the lattice points 

separated by a distance a apart when the kink velocity is approximately 

constant which is Q = 2nX/a = KX, which we call the kink frequency, and 

its harmonics Qj = 2irjX/a = KjX. Which can also be seen by setting 

X(t) ~ Xt in sin27r/'X(£)/a which becomes sin 27r(KjX)t = sin27rlX,-£. Note 

this frequency is determined only by the periodicity of the lattice and does 

not depend on the strength of interaction as Upn does. The Kj are the 



vectors of the reciprocal lattice. For radiation of frequency f2 the kink has to 

be essentially free in the sense that the kinetic energy is appreciably greater 

p 2 
than the potential energy i.e., -^- ^> Vpn- In order to see the role played 

by Q in kink radiation processes consider the solution for 

rt 
(j) n {t) = J2*l>k(n) exp[-iiu(k)t]J2C jk / dt' expi[cu(k)t' -2njX(t')/a] (2.7) 

3 ° 

where 

n 

Now if we can treat the kink as a free particle for times t small compared with 
the radiation relaxation time but long with a suitably short interaction time 
we can replace X(t) by Xt in Eq. (2.7) and obtain the usual delta function 
of perturbation theory exactly as in Fermi's golden rule. In the large system 
the phonon levels become continuous and the interaction time is Uq 1 . When 
we let t — > oo in Eq. (2.7) we obtain 

0(n, t) = I dktpk(n) exp[— iu(k)t] ^^Cjk 5[u(k) — 2-irjX/a] 

3 

= ^^kU)(n)exp[-i2irjX/a]tC jk (j)u[k(j)]/k{j) (2.8) 

3 

where k(j) is the solution for k of the radiation resonance condition 
Qj = u[k(j)}, and where for convenience we have taken units where 7r// = 1. 
Thus for a given j, say j = 1, our solution, (f>(n,t) is a single normal mode 
V'fe(i)( 71 ') with a single frequency 2ixX/a which is the phonon frequency ra- 
diated by the kink which is resonant with the kink frequency KjX. The 
general solution is a linear combination of those phonons in the phonon band 
which are resonant with the kink frequency and its harmonics KjX. The 
damped driven small system does not have a delta function radiation reso- 
nance condition but a Lorentzian centered on the resonance condition with 
a finite linewidth. We discuss the effect of the finite linewidth on Sees. Ill 
and IV. The apparent divergence of Eq.(2.8) at k = is spurious and is not 
a breakdown of perturbation theory. It is a breakdown of the assumption 
that X(t) is a constant which is true everywhere except at the band edge 
where the time dependence of X(t) has to be considered within perturbation 
theory. For a full discussion of the treatment of X(t) in the large system 



see [4]. The spurious divergence does not arise in the small discrete system 
because of the external damping T<p n . 

For each phonon mode, ui(k), there will be a set of velocities X such 
that Qj = u(k). The ratio of successive kink velocities resonant with a 
given u(k) is (Xj + \/Xj) = (j/j + 1). The number of harmonics Qj and the 
number of kink velocities that correspond to a given phonon mode, u(k) is 
determined by the discreteness parameter (tt/Iq). As the system becomes 
more discrete i.e., as Iq decreases the more Cj will contribute to Eq. (2.3) 
for Vpn resulting in more harmonics Qj contributing to the radiation by 
the kink. However as the system becomes more discrete the PN potential 
well will become deeper and the kink becomes trapped sooner, the range of 
velocities for which ^f- ^> Vpn decrease leading to an upper limit on the j 
values such that the "free condition" is valid. When fij lies in the phonon 
band the kink loses energy by radiation of phonons. The unperturbed energy 
of the kink plus phonons is conserved. When the velocity of the kink decays 
to Xj = (2irj)~ 1 uo the kink abruptly stops radiating at the frequency IX,- 
because there are no phonon states u(k) < ujq. This sharp decrease in the 
rate of loss of energy by the kink is called a knee in Ref . [3] . The remaining 
radiation comes form higher harmonics Qj> = 2nj f X /a where j' > j which 
are still in the band and still radiating. When the velocity of the kink decays 
further to X,j + \ = [2n(j + l)] -1 ^ the kink stops radiating the j+1 harmonic 
causing the next knee in the radiation rate. The pattern is repeated until 
the kink loses enough energy to be trapped in a PN well where it will radiate 
the harmonics of ojpn at a much lower rate. 

Peyrard and Kruskal [3] simulated and analyzed the discrete SG for a 
large system and found for a variety of parameter values (tt/Iq = .9, 1, and 
1.05) and initial values of the velocity that there were one, two or three 
knees. They obtained numerical values for the velocity knees by a graphical 
solution of the discrete SG dispersion law which agreed with their simulations 
to better than one percent. Their numerical values for the three knees are 
given by our radiation resonance condition for j = 1,2 and 3 to within one 
percent. The velocity of the jth knee is given by Xj = (2tij)^ 1 ujo as was 
pointed out in [2]. 

There is one continuum SG kink radiation problem that gives further 
insight into the general kink radiation problem and that is the model of the 
continuum SG with a spatially periodic perturbation potential which was 



analyzed in Ref. [4]. The perturbed SG equation is 

<j> - <p" + (tt//o) 2 (1 + e sin k'x) sin = 

where e is a small dimensionless parameter and k' is the wave number of the 
perturbation. The competing length scales are Iq 1 and a i.e., the size of the 
kink and the lattice spacing. When we go to a CV description [4] by making 
the ansatz = <t[(ti/1 )(x — X(t))] + x[( n /^o)( x — X(t)),t] where a is the 
kink solution of Eq. (2.8) when e = and to lowest order in e we obtain the 
following coupled equations of motion 



and 
where 



X ~ x" + (tt/^o) 2 cos ax = -e(ir/l ) 2 sin k'X 



X + u 2 sin k'X = 



uj 2 = (e/l )(k'l /2) 2 [sm(k%/2)]-\ 



The equation for X is the pendulum equation. In [4] the radiation and 
X(t) are calculated analytically. The relevance of this example for the present 
discrete case is that the X equation is of the same form as Eq. (2.3) with 
C\ = iJ 1 and Cj = for j > 2. Consequently only one velocity knee occurs 
i.e., when = X = (uo/kf) which is the same as our radiation resonance result 
for K\ = k'. The kink radiates phonons at frequency ui(k) and wave number 
k such that k'X = u(k). The example of the continuum SG kink in the 
(esink'x)(f) perturbation makes clear that the kink frequency which is reso- 
nant with the phonon frequency u(k) when X(t) ~ Xt is determined by the 
spatial periodicity of the perturbing potential which is K~ l in the discrete 
case and (k')~ l in the perturbed continuum case. The resonance condition 
arises in the same way in both cases i.e., 5[u(k) — k'X] and S[u(k) — KjX] 
after setting X(t) ~ Xt. We next consider the kink in small finite sys- 
tems which satisfy periodic boundary conditions as in the experiments with 
Josephson junction rings, [1], where cases of 4 and 8 degrees of freedom with 
1,2,3 or 4 kinks are treated. In the small discrete periodic systems the kink 
interacts with the radiation it has spontaneously emitted. However, the in- 
duced absorption and emission processes are much smaller than the dominant 
spontaneous emission process and don't play a significant role. They have 
the same radiation resonance condition as the spontaneous emission process. 



One difference appears in the small periodic system is the phenomena of 
hysteresis which we discuss in Sec. IV. 

In the Josephson junction rings experiments the kink radiation interac- 
tions are perturbative which implies the phonon eigenmodes and eigenvalues 
are those of the linear phonon equations in the presence of the kink. Note 
the number of modes and the spectra u(k) are unaffected by the presence of 
the kink. Since many discussions of the experiments often have an incorrect 
number of modes and frequencies of the modes it is worth while to provide 
a brief elementary review [5]. The phonon modes in a lattice of spacing a 
with N particles satisfying periodic boundary conditions have displacements 
u(sa) where s is an integer which satisfy u[sa] = u[sa + L] = u[(s + N)a] 
where L = Na. The solution for u is u[sa,t] = M [o]e*' sA:a " a; ^- )t l where 
k = 0, ±2tt/L, ±Aix/L ■ ■ ■ Nir/L and u(k) = [1 + 4 sin 2 (ka/2)] 1 / 2 and where 
for convenience we use units where (vr// ) = 1. 

The number of phonon modes is equal to the number of particles N which 
is also the number of k values i.e., k — 0, ± • • • ± (N — 2)ir/L, Ntt/L where 
the modes k = and Ntt/L are non degenerate and the remaining modes are 
two fold degenerate with one mode traveling in the clockwise direction and 
one in the counterclockwise direction with the same frequency uu(k) = u(—k). 
The k = mode corresponds to all N particles oscillating in phase at the 
frequency c^o- The k = tt corresponding to A = 2, where each particle 
oscillates 180 degrees out of phase with its nearest neighbor at a frequency 
of (uq + 4) 1 / 2 . The number of different phonon frequencies u(k) is [(N — 
2)/2 + 2] = (JV/2) + 1. For a typical case of iV = 8 there are five different 
frequencies u;(0), a;(7r/4), lj(tt/2), a;(37r/4), and uj(tt) where we set a — 1. 
The crucial point is that in the radiation resonance condition, i^j X = u(k), 
that the Kj are vectors of the reciprocal lattice and are completely different 
than the vectors k that determine the phonon modes of the lattice. As an 
illustration consider the case N = 8, j = 1 and M — 1. For each of the five 
different values u>(k) there will be one velocity of the kink which is resonant 
namely, X(k) = a(27i)~ l uj(k) for a total of five different velocities. The next 
question that arises is given the number of phonon frequencies is fixed by 
N i.e., (N/2) + 1, how many different kink velocities can be resonant with 
the five modes. The answer is determined by the number of j values that 
contribute to Vpn without the kink being trapped. For each j value there 
are [(N/2) + 1] velocities so the answer is j[(N/2) + 1] different values of X 
and thus the same number of voltages would be possible. The number of j 
values is determined by the discreteness parameter tt/1 . 

10 



There is a theorem [6] which relates the total number of degrees of freedom 
to the number of kinks and phonons which states the total number of degrees 
of freedom N must satisfy N = Nk + / p(uj)duj where Nk is the number of 
kinks and / p(u)duu is the effective number of phonon degrees of freedom 
with p(u) the phonon density of states. This theorem does not reduce the 
number of eigenfrequencies, u a , but their effective weight i.e., in our case 
p{u>) = E a Ca^( w ~ w a ) thus we have N = N K + J^C a . The states u a are 
the same as in the absence of the kink but the effective number of phonon 
states is J2C a = N — Nk reduced from the N it would be in the absence of 
kinks. 

In this section we have seen the radiation resonance condition at line 
center, KjX = u(k) is the same in the large discrete periodic system. The 
number of j values i.e. harmonics of 2nX /a is determined by the boundary 
conditions and the dispersion law for u(k) is the same in both cases. The 
small periodic system requires external damping to justify radiation pertur- 
bation theory while in the infinite system the quasi continous density of states 
is sufficient. We defer to Sec. IV a discussion of hysteresis which occurs in 
the small periodic system and for which there is no analogue in the large 
system. 

III. ANALYSIS OF JOSEPHSON JUNCTION RING 
EXPERIMENTS 

In a recent paper Watanabe et al. [1] presented an analyzed experimental 
measurements and simulation results for the current- voltage (IV) character- 
istics of rings of 4 and 8 underdamped Josephson junctions using a model of 
the damped driven SG equation with periodic boundary conditions. The SG 
equation in their notation is $ n +r<l> n +sin $ n — A 2 A 2 $ n — I /Ic where V is the 
damping and I/I c is the damped driven form of Eq. (2.2) the dimensionless 
parameter A 2 is associated with the A 2 term while in Eq. (2.2) the param- 
eter is associated with the nonlinear sin0 n term. They found two types of 
traveling waves low velocity kinks and high velocity whirling modes. In this 
paper we are concerned only with the low velocity kink modes. Their analy- 
sis of the behavior of the low velocity kinks in their words "consists of rough 
physical arguments that are far from rigorous." In this section we show the 
CV radiation perturbation theory in Sec. II correctly explains the number of 
peaks and their position in the I-V graphs very accurately. In particular the 
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experimental peaks are determined by the phonon normal modes, the kink 
frequencies, 2njX /a, the radiation resonance, condition, 2ttjX /a = u(k), as 
presented in Sec. II. 

In order to explain the I-V diagrams we first review the voltage phase 
relationship as given in Eq. (2.16) of [1] as the Josephson relationship Vj = 
Vo&j/A where for a kink $., = &jX which shows for a kink traveling with 
a constant velocity that Vj for each j is ~ X. For a phonon normal mode 
all the <t>j ~ u{k) and thus we have Vj ~ uj{k). For a kink moving at 
constant velocity X = a/T, where T is the time to travel the lattice spacing 
a, we have the kink frequency u = 27rX/a in the dimensionless units of 
Eq.(2.2). Thus the space and time average of Vj/Vq = V/V for a constant 
velocity kink is V/Vq = 2nX/a or V/Vq = u/A in the dimensionless units of 
[2.2]. In the case of a phonon frequency we have V/Vq = u(k)/ A = [A~ 2 + 
4sin 2 (/ca/2)] 1/ ' 2 = [u; 2 + 4sin 2 (/ca/2)] 1/ ' 2 where in the last equality we used the 
relationship between the two dimensionless set of variables of [2.2] and [3.1] 
where (tt/1 ) 2 = A~ 2 . However because of the radiation resonance condition 
the frequency of the phonon and thus V/Vq = 2irX /a = [c^Q+4sin 2 (A;a/2)] 1 / 2 
which is true independent of space time averaging. 

Thus the line center of each phonon line, u(k), determines a unique value 
of the voltage V. Since X is a nonlinear function of / we have an I(V) 
for each phonon line. We show how the finite phonon linewidth in a small 
discrete periodic systems leads to the full I(V) lines and hysteresis later in 
this and the next section. We also point out that the preceding conclusion for 
V/Vq = 27rX/a = uj(k) is valid also for the harmonics of the kink frequency 
i.e. V/Vq = 2njX /a = ui{k) for j an integer as long as for those j values the 
kink is untrapped as discussed in Sec. II. 

In order to avoid all effects of space time averaging I compare the ra- 
tio of the experimental voltages to the ratio of the corresponding phonon 
frequencies and thus to the ratio of the kink velocities. I now show the 
CV kink-radiation theory of [2,3] determines the number and position of the 
voltage lines in the I-V diagrams for all the experiments and simulations 
including new results for the multi kink cases. We consider the individual 
cases one at a time. The first case we consider is a single kink, M — 1, 
N = 8 junctions, discreteness parameter A = 1.48 and T = 0.13 in fig.10 
of [1]. In fig. 10 there are five peaks in the I-V diagram determined by the 
five minima of the mi/i ) curve which corresponds to the correct five phonon 
frequencies for N = 8 junctions with periodic boundary conditions. (On Fig. 
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10 a sixth point is indicated to be a minimum but does not appear to be 

from the figure.) We next show that the measured voltage ratios agree with 

the corresponding, phonon frequency ratios to about 2.5%. The ratio we 

calculate are 

uj(k + l)/u{k) - V(k + 1)/V(k) 



ou(k + l)/u(k) (3 ' 2) 

where the u){k) are given by Eq.(2.5) for k = 0, ±7r/4, ±7r/2, ±37r/4, n and 
the voltage ratios are obtained from the experimental values given in Fig. 
10. The ratios of the succesive calculated phonon frequencies are 

^M = 1.513, ^ = 1.54, ^/i>=1.26, -#^ = 1.07 
w(0) cj(tt/4) ' u(ir/2) ' cu(3tt/4) 

while the corresponding experimental values of the four voltages are 
V{1) - 1 47 V{2) - 1 54 V(3) - 1 29 V(4) - 1 107 

no)- 1 - 47 ' W)~ ' v(2) _ ' nay - 

where V^(0) is the voltage corresponding to w(/c) for A; = and the minima 
of j/j/ f c ? , ^(1) is the second minima of i-//'? corresponding to cj(7t/4) 
etc. When we substitute the phonon frequencies a; (A; + \)/uj{k) and voltages 
V(/c + 1)/V(k) in Eq.(3.2) we obtained for the relative percentages 2.8%, 
2.6% and 3% for an average 2.5% which demonstrates impressive agreement 
with the resonance condition 2irX/a = u(k). Our results confirm that when 
V(k)/Vo is plotted versus the five k values you obtain the function [cj 2 + 
4sin 2 (/ca/2)] 1//2 evaluated at the corresponding k value as you should. The 
discreteness parameter in Fig. 10, A = 1.48, is not sufficiently small enough to 
cause the discreteness required to see the j '• = 2 harmonics in 27rjX/a = u(k). 
However when Watanabe et al simulated the I — V curve with M — 1, N — 8, 
T = .02 and a smaller and thus more discrete A 2 = 1 there were indications 
of two j = 2 lines in their Fig. 11. The voltage of the line they labeled 
m = 15 is one half of the voltage of the line m = 7 and the voltage of their 
m = 9 line is one half of the m = 6 line. The agreement of the ratios is 
less than three percent. If the two lines m = 15 and m = 9 are j = 2 lines 
than there are only three j '• = 1 lines instead of the 5 j = 1 lines remaining. 
When we assign the three remaining lines as follows o;(0) ^ line 8, uj(tt/A) ^ 
line 7 and uj(3tt/A) ^ line 6 their ratios compare with the simulation ratios 
for Vg, V7, and Vg and voltages in Eq. (3.2) to better than three percent. 
Furthermore the frequencies a;(37r/4) and uj(tt) are within six percent of each 
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other and Fig. 11 is such that the I and V values are not carried out far 
enough to be able to see the full uj(3tt/4) and uj(tt) lines. Assuming that 
this is correct then we have 4 j = 1 lines and 2 j = 2 lines that agree well 
with the radiation resonance condition 2irjX = uj(k). There is one problem 
and that is the missing j = 1 line, the uj(tt/2) line does not appear on the 
numerically calculated I-V curve in Fig. 11. Its absence is unexplained and 
consequently all I can conclude is that we are probably seeing 4 j = 1 and 2 
j = 2 lines. 

In addition to performing experiments on iV = 8 junction systems Watan- 
abe et al did a thorough study on an N = 4 system with one kink, A 2 = 1 
and r = 1. The allowed k values are k = 0, ±7r/2 and it with frequencies 
u(0) = 1, cu(tt/2) = 1.73, and u(ir) = 2.24. For N = 4 there should be and 
are three frequencies in their Fig. 9. When we take the measured voltages in 
Fig. 9 and form the ratios in Eq.(3.2) we obtain the first ratio agrees to within 
two percent and the second ratio is essentially exact so the average agreement 
with the phonon resonance condition is one percent. The authors in [1] pro- 
vide further insight into the above N — 4 system with the same parameters 
in their Fig. (8) by plotting (f>j(t) vs. t for each of the four junctions. From 
the figure you can see that the (f)j (t) is a superposition of a kink and a phonon 
with wavelength A = 4a, where 4a is the size of the system, which means the 
k of the phonon is tt/2 because ka = tt/2 =>- A = 4a. Consequently the veloc- 
ity of the kink is X/a = (2tt)- 1 uj(tt/2) = (2tt)- 1 (3) 1/2 . We can represent the 
kink and phonon configuration as ^x ^^^- as it appears in Fig.8 where x 
is the center of the kink and the arrows indicate the phase of the k = tt/2 
phonon mode at the four junctions, and there is a change of of 2tt over 
the four junctions. Thus Fig. 8 is a beautiful conformation of the fact that 
the kink and the phonon have exactly the same frequency 2irX/a = uj{tt/2) 
as they travel around the ring. Next we consider the multiple kink case and 
find a completely new effect in which the spacing of kinks in the ring allow 
only certain k value for the resonant phonons which we can explain by the 
CV radiation theory. The case of multiple kinks treated in Fig. 4 for the 
parameters N = 8, A 2 = 4, r=.17, and for M=l,2,3 and 4 kinks. The value 
A 2 =4 (which is only weakly discrete) corresponds to a kink size, 47ra which 
is larger than 8a, the size of the system, which means all the kinks overlap 
each other. In the continuum the SG is integrable and thus the kinks do not 
interact with each other. In the present weakly discrete case the discreteness 
is effectively a perturbation which destroys the integrability and each kink 
interacts with all other kinks weakly because they all overlap each other. It 
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would be a difficult problem to work out the solution of the equations of 
motion for TV interacting radiating kinks traveling around the ring. However 
as we see below we can infer the kinks travel around the ring keeping equal 
spaces between each other. What is important is that the phonon modes 
are independent of the number of kinks and are the same as the modes in 
the absence of the kinks which is a consequence of the fact that kink-phonon 
interaction in all cases in this paper can be treated perturbatively. For the 
present case, A 2 =4, the kink phonon interaction is sufficiently weak that only 
the lowest harmonic, j — 1, is observed. For N = 8 as before there after five 
frequencies that correspond to the k values k — 0, ±7r/4a, ±7r/2a, ±37r/4a, 
and ±7i /a. 

We have as before (V/V ) = (u/A) = [A~ 2 + 4sin 2 (A;a/2)] 1 / 2 = [^ + 
4sin 2 (/ca/2)] 1 ' 2 = (2%MX /auJa). The fourth equality is the radiation reso- 
nance condition for M kinks. The only changes for M > 2 from the M = 1 
case is that the phase change for M kinks is equal to 2ttM, (2tt) for each 
kink and the corresponding dimensionless frequency is 27rMX/auo, where 
all kinks have the same velocity. The total rate of change of phase for a 
given X is M times the rate of change of one kink, i.e., 2ttX /oujq. 

In the experiment reported in Fig. 4 of [1] there are a total of four lines 
one for each of the four values of M i.e., M = 1,2,3,4. When we make 
the assignment V(M)/V = [cu 2 + 4sin 2 (A; M a/2)] 1/2 = 27iMX M /uj a where 
k\ = 7r/4a, k 2 = vr/2a, k% = 37i~/4a and k^ = ir/2a we find the experimental 
values of V(M + 1)/V(M) agree with the calculated frequency ratios oj{k,M + 
1)/u;(/cm) to within two percent. Consequently the frequency assignments 
uj{M) are clearly correct. The question arises why is there only one frequency 
for each value of M and why is that the frequency corresponding to the 
assigned values k 2 = (7i/2a),k 3 = (Sir /4a), and k A = (n/a). We find the 
answer by finding the wavelength Am where ku = 2n/\ M . For M = 2 
A2 = 4a, for M = 3 A3 = 8a/3, and for M = 4 Am = 2a. Thus we find 
the wavelength of the resonant phonon for each M resonance is exactly the 
spacing between each other of the M kinks assuming the kinks are evenly 
spaced i.e., for M = 2 kinks the spacing between the kinks is 4a, M = 3 the 
spacing is 87r/3 and M = 4 the spacing is 2a in agreement with the kM values. 
(Recall the total length of the ring is 8a) The explanation of the results are 
that the kinks are the source of radiation in the linear inhomogenous wave 
equation for the phonons. We have already seen that the frequencies in the 
source are the phonon frequencies 2-nMXja for each M but when we have 
M > 2 we also have a spatial dependence in the source namely the equal 
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spacing between the kinks which determine the k values of the phonons. That 
is the spatial properties of the phonons radiated by the kinks are determined 
by the spatial dependence of the source which is given by the arrangement 
of the M kinks rotating around the ring. The only assumption is that for 
M > 2 the kinks moving around the ring are evenly spaced. One would 
have to solve the full M interacting kinks and radiation problem to prove 
the assumption. Alternatively one could reasonably conclude that in the 
solution the kinks were evenly spaced by the agreement of the experimental 
results with the assumption of equal spacing. 

Up to point we have been using the radiation resonance condition to 
determine the voltage at the line center of the line on the I — V diagram. In 
the infinite discrete SG equation the phonon band is a continuum and the 
radiation resonance condition 27rjX/a = ujq picks out a definite value of X 
because the phonon line is sharp. In the small system case we have a discrete, 
damped and driven SG, Eq.(3.1), so that the valid perturbation theory leads 
to a Lorentzian line shape for the phonon line which in turn generates the step 
in the I — V curves. There are two causes for the linewidth in small discrete 
systems. First the damping term Tipj gives a contribution to the linewidth 
proportional to T and contributes to each mode in the same way. The second 
mechanism is the radiation linewidth due to spontaneous emission, r~ 1 (k), 
which is different for each k thus causing the linewidth of u(k) to depend on 
k. In addition, r~ 1 (/c) is a nonlinear function of X. As a result, the delta 
function that appears in the infinite system case 5[2njX/a — oo(k)] becomes 
a Lorentzian in the damped driven small system, case 

[{2irjX/a - Lu(k)) 2 + (T + r' 1 (k)) 2 ]- 1 (T + r" 1 ^)) (3.3) 

The voltage V is proportional to X which is a nonlinear function of / 
through r(k). Consequently in the small system V is a nonlinear function 
of / which leads to the nonlinear I — V curves in the damped driven small 
system. The curves in the I—V diagrams for different ui(k) thus have different 
widths and the lines from adjacent u(k) overlap each other leading to a step 
structure and hysteresis. We explain the step structure and hysteresis in 
more detail in the next section. 
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IV. DISCUSSION AND CONCLUSIONS 

In summary we have shown when the kink-phonon interaction can be treated 
perturbatively that we can treat the discrete SG or (any discrete non linear 
Klein-Gordon equation) as a "point particle," X(t) interacting with a gas 
of non interacting phonons where the effective frequency of the kink is Q = 
2-nX ja = KX and its harmonics Qj = jQ where j is an integer. The 
resonance conditions are Qj = u(k). In the large system the phonon band 
is effectively a continuum and the resonance condition appears as a delta 
function S[Qj — u(k)]. In the small system we have a Lorentzian with the 
same resonance condition and a linewidth (r + r~ 1 (A;)). For the perturbation 
theory to be valid the system must be discrete enough for a moving kink to 
produce radiation but not so discrete that the kink gets trapped in a PN well 
or that nonlinear phonon-phonon interactions become important. 

In particular we have explained the experiments and simulations of a small 
discrete periodic dynamically arrays of Josephson junctions in [1]. Watan- 
abe et. al [1] ended their paper by listing three open problems requiring 
further research in the low voltage region of small discrete Josephson junc- 
tion arrays. In the present paper we have solved all three problems. Their 
first problem requiring solution was can one find the correct resonance con- 
dition and derive it. The correct answer which we obtained in this paper 
is, 2irjX/a = Qj = u(k), i.e., the kink frequency is equal to the phonon 
frequency which is the usual radiation resonance condition that comes from 
the first Born approximation. 

Their second problem to be answered was to find a precise characteriza- 
tion of the steps in the low voltage region of the I — V diagram. The steps in 
the I — V diagrams are a direct consequence of the kink phonon interactions. 
The separation of the phonon states in the small Josephson junction rings 
is comparable or greater than their width T + r~ 1 {k). Since r depends on 
k the width of the steps are different for different k. Consider the driver / 
driving the kink at a velocity X such that the kink frequency 2i\X ja is in 
the region roughly between the phonon frequencies uj(k — 1) and u(k). Then 
practically all the energy from the driver goes into increasing the velocity 
from the kink. As we increase I then X (which is proportional to the voltage 
V) increases until eventually 27iX/a comes into the linewidth of u(k) and as 
the kink frequency approaches u(k) the kink begins to radiate. Thus / has 
to increase rather sharply to provide energy to the kink to compensate for 
the energy loss of the kink due to its radiation of phonons and to prevent 
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X from decreasing. Consequently there is a sharp increase in I for approxi- 
mately the same V i.e., X. As we increase / further X starts to increase thus 
it leaves the center of the u(k) line where the radiation from the u(k) line 
decreases to zero. Now as / increases nearly all the energy from I goes into 
again increasing X until the kink frequency enters the wing of the u(k + l) 
line and the whole process repeats itself. Thus we have a precise qualitative 
and quantitative characterization of the steps. 

In the previous paragraph for convience I discussed the steps as if the 
linewidths of u(k) and u(k + 1) had a small overlap in the wings of the lines. 
However the wings of the lines can overlap a nontrivial amount and the ex- 
planation of the formation of the steps remains valid but a new phenomena 
appears namely hysteresis. We can understand the phenomena of hysteresis 
by the following argument. In each step u(k) the driver is a different non- 
linear function of voltage because r(k) depends on both k and X. Or stated 
differently the non linear function I(V) is a different function for each allowed 
k. When the linewidths of u(k) and u(k + l) overlap in their respective wings 
there are two stable states with the same X, i.e., voltage for a given / corre- 
sponding to the r{k) of u(k) and r(k + 1) of cu(k + 1). Starting at the line 
center of u(k) as / is increased X increases and the u linewidth becomes less 
stable and the u(k + 1) linewidth becomes more stable until the stability of 
the u(k) vanishes and the system jumps to the u(k + 1) line. Instead if you 
start at the line center of u(k + 1) with linewidth r(k + 1) and decrease I 
then X decreases and eventually the system jumps discontinuously from the 
u(k + 1) line to the u(k) line at a different value of / (than the increasing / 
case) because r(k + 1) depends differently on X then r{k) does. 

Their third problem requiring solution is "can the prediction of the high 
wave number resonances be experimentally observed?" As we discussed thor- 
oughly in Sec. II, the resonance 27ijX/a for j'=2,3.. are observable if the 
system is discrete enough i.e., if (tt/Iq) is sufficiently large or equivalently 
when A is sufficiently small. The more discrete the system the more j values 
that are observed. There is a maximum j determined by the discreteness 
that occurs when the kink becomes trapped in a PN well. 

We conclude with the observation that the explicit solution for the steady 
state radiation intensity and the explicit solution for r(k, X) can be obtained 
by solving for the steady state solution of the linear Eq. (2.4) and calculating 
the radiation relaxation time, r, in first order perturbation theory by the 
methods in [4] and [7]. 
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